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0.1 INTRODUCTION 

Since the early works [1-4] on the so-called nondiffracting waves (called also Localized 
Waves), a great deal of results has been published on this important subject, from both 
the theoretical and the experimental point of view. Initially, the theory was developed tak- 
ing into account only free space; however, in recent years, it has been extended for more 
complex media exhibiting effects such as dispersion[5-7], nonlinearity[8], anisotropy[9] 
and losses[10]. Such extensions have been carried out along with the development of ef- 
ficient methods for obtaining nondiffracting beams and pulses in the subluminal, luminal 
and superluminal regimes[l 1-18]. This paper (mainly a review) addresses some theoreti- 
cal methods related to nondiffracting solutions of the linear wave equation in unbounded 
homogeneous media, as well as to some interesting applications of such waves. 

The usual cylindrical coordinates (p, cf), z) will be used here. In these coordinates the 
linear wave equation is written as 

p dp \ dp ) p 2 d4> 2 dz 2 c 2 dt 2 

In section II we analyze the general structure of the Localized Waves, develop the so 
called Generalized Bidirectional Decomposition, and use it to obtain several luminal and 
superluminal (especially X-shaped) nondiffracting wave solutions of eq. d0.lt . 

In section III we develop a space-time focusing method by a continuous superposition 
of X-Shaped pulses of different velocities. 

Section IV addresses the properties of chirped optical X-Shaped pulses propagating in 
material media without boundaries. 

Finally, in Section V, we show how a suitable superposition of Bessel beams can be 
used to obtain stationary localized wave fields with a static envelope and a high transverse 
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localization, and whose longitudinal intensity pattern can assume any desired shape within 
a chosen interval < z < L of the propagation axis. 

0.2 SPECTRAL STRUCTURE OF THE LOCALIZED WAVES AND THE 
GENERALIZED BIDIRECTIONAL DECOMPOSITION 

An effective way to understand the concept of the (ideal) nondiffracting waves is furnishing 
a precise mathematical definition of these solutions, so to extract the necessary spectral 
structure from them. 

Intuitively, an ideal nondiffracting wave (beam or pulse) can be defined as a wave capable 
of maintaining indefinitely its spatial form (apart from local variations) while propagating. 

We can express this intuitive property saying that a localized wave has to possess the 
property[12,13] 

A 7 

V(p,(t>,z,t) = *(p,0,z + Az o ,i+-^) (0.2) 
where Az is a certain length and V is the pulse propagation speed that here can assume 
any value: < V < oo. 

Using a Fourier Bessel expansion, we can express a function ^>(p, <j>,z,t) as 



V(p,<l),z,t) = 



n=— oo 



oc 

'/ dk p I dk z I d^k p A n {k p ,k z ^)J n {k p p)^ z e- iut ^ 
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(0.3) 

Using the translation property of the Fourier transforms T[f(x + a)] = exp(ika)T[f (x)], 
we have that A n (k p , k z ,uj) and exp[i(k z AzQ — u>Azo/V)]A n (k p , k z ,u>) are the Fourier 
Bessel transforms of the l.h.s and r.h.s. functions in eq. d0.2l l. And from this same equation 
we can get[12,13] the fundamental constraint linking the angular frequency uo and the 
longitudinal wavenumber k z : 

ui = Vk z +2rrm-— (0.4) 
Az 

with to an integer. Obviously, this constraint can be satisfied through the spectral functions 

A' n (k p ,k z ,uj). 

Now, let us explicitly mention that constraint (10. 4t does not imply any breakdown of the 
wave equation validity. In fact, when inserting expression (10.3b in the wave equation (??), 
one gets that 

LO 2 



= K + K (0.5) 

So, to obtain a solution of the wave equation from (10.3b . the spectrum A n (k p , k z ,u>) 



must have the form 



A n (k p ,k z ,u) = A n (k z ,Lu)S[k 2 p -[^-ki)\ (0.6) 
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where S(.) is the Dirac delta function. With this we can write a solution of the wave equation 
as 



n— — oo 



dto I dk z A n (k z ,Lu)J n pyj-g - kl e* 2 e- mi e m * 

(0.7) 

where we have considered positive angular frequencies only. 

Equation ( 10.7b is a superposition of Bessel beams and it is understood that the integrations 
in the (to, k z ) plane are confined to the region < w < oo and — oj/c < k z < lo/c 

Now, to obtain an ideal nondiffracting wave, the spectra A n (k z , oj) must obey the fun- 
damental constraint ( 10.4) . and so we write 

oo 

A n (k z ,u) = J2 S nm (u,)5(Lu~(Vk z +b m )) (0.8) 

m— — oo 

where b m are constants representing the terms 2mirV/ Azo in eq.( 10.4l i. and S nm (u>) are 
arbitrary frequency spectra. 

With eq.( 10.8t into eq. (IQ.71 >. we get a general integral form of an ideal nondiffracting wave 
defined by eq. (10.2b : 

oo oo 

*(p,0,«,t) = 1Pn m (P,<l>,Z,t) (0.9) 

n— — oo m— — oo 

with 



ipnm(p,<t>,Z,t) = e ~ lbmZ/V / dw S nm (u) 



(0.10) 

where Lo m - m and w max depend on the values of V: 

• For subluminal (V < c) localized waves: b m > 0, (w m i n ) m = cb m /(c + V) and 

(Wraax)m = cb m /(c - V). 

• For luminal (V = c) localized waves: b m > 0, (u> m m)m = b m /2wd (w m ax)m = oo- 

• For superluminal (V > c) localized waves: b m > 0, (w m in)m — cb m / (c + V) and 
Kax)m = oo. Or b m < 0, (w min ) m = cb rn j (c - V) and (w max )m = oo. 

It is important to notice that each ip n m{p, 4>, z, t) in the superposition (10.9b is a truly 
nondiffracting wave (beam or pulse) and the superposition of them, (10.91 1. is just the most 
general form to represent a nondiffracting wave defined by eq. (10.2l >. Due to this fact, the 
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search for methods capable of providing analytical solutions for ip nm (p, <p, z, t), eq. d0.10b . 
becomes an important task. 

Let us remember that equation (10.101) is also a Bessel beam superposition, but with the 
constraint (10.41) between their angular frequencies and longitudinal wavenumbers. 

In spite of the fact that the expression (10.101) represents ideal nondiffracting waves, it is 
difficult to obtain closed analytical solutions from it. Due to this, we are going to develop a 
method capable of overcoming this limitation, providing several interesting localized wave 
solutions (luminal and superluminal) of arbitrary frequencies, including some solutions 
endowed with finite energy. 

0.2.1 The Generalized Bidirectional Decomposition 

For reasons that will be clear soon, instead of dealing with the integral expression (10.9I) . our 
starting point is the general expression (I0.7l i. 

Here, for simplicity we will restrict ourselves to axially symmetric solutions, assuming 
the spectral functions so that 

A n (k z ,oj) = S n0 A(k z ,u;) (0.11) 

where S n o is the Kronecker delta. 

In this way, we get the following general solution (considering positive angular frequen- 
cies only) which describes axially symmetric waves: 



duj 



dk z A(k z ,uj)Jo p\ -it 



— U)f c 




(0.12) 



As we have seen, ideal nondiffracting waves can be obtained since the spectrum A(k z ,u>) 
satisfies the linear relationship ( 10.41 ). In this way, it is natural to adopt new spectral param- 
eters in the place of (u>, k z ) that make easier to implement that constraint[12,13]. 

With this in mind, we choose the new spectral parameters (a, (3) BY 



2V 



(w - Vk z ) 



(0.13) 



Let us consider here only luminal (V = c) and superluminal (V > c) nondiffracting 
pulses. 

With the change of variables (10. 13b in the integral solution d0.12l l. and considering ( V > 
c), the integration limits on a and (3 have to satisfy the three inequalities 



< a + (3 < oo 

c+V 
c + V n 



(0.14) 



a > 
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Let us suppose both a and /3 to be positive [a, (3 > 0]. The first inequality in (10. 14b is then 
satisfied; while the coefficients (c— V) / (c+V) and (c+V) / (c— V) entering relations (10. 1 4b 
are both negatives (since V > c). As a consequence, the other two inequalities in (10.141) 
result to be automatically satisfied. In other words, the integration limits in < a < oo 
and < (3 < oo are contained in the limits (10.14b and are therefore acceptable. Indeed, 
they constitute a rather suitable choice for facilitating all the subsequent integrations. 

Therefore, instead of eq.( 10.12b . we shall consider the (more easily integrable) Bessel 
beam superposition in the new variables [with V > c] 



jf daj^ d[3A(a, /?) Jo (p\j{^- l) (a 2 + /9 2 ) + 2 (Yl + lj ap\ 



(0.15) 

where we have defined 



C = z-Vt; r) = z + Vt (0.16) 

The present procedure is a generalization of the so-called "bidirectional decomposition" 
technique[l 1], which was devised in the past for V = c. 

From the new spectral parameters defined in transformation ( 10.731 ). it is easy to see that 
the constraint ( 10.41 ). i.e. lu — Vk z + b, is implemented just by making 



A{k Zl uo)^ A(a,(3) = S(a)6{(3-(3 ) (0.17) 

with fto = b/2V. The delta function 5((3 — 0o) in the spectrum (10.17l i means that we are 
integrating Bessel beams along the continuous line u> = Vk z + 2Vfio and, in this way, the 
function S(a) will give the frequency dependence of the spectrum: S(a) — > S(uo/V — (3o), 
This method is a natural way of obtaining pulses with field concentration on p — and 

Now, it is important to stress[13] that, when /3q > in (10.17b . the superposition (10.15b 
has contributions from both backward and forward Bessel beams in the frequency intervals 
VPo <ui< 2VPo (where k z < 0) and 2V(3q < ui < oo (where k z > 0), respectively. 

Nevertheless, we can obtain physical solutions when making the contribution of the 
backwards components negligible, by choosing suitable weight functions S(a). 

It is also important to notice that we use the new spectral parameters a and [3 just to obtain 
(closed-form) analytical localized wave solutions, IT being that the spectral characteristics 
of these new solutions can be brought into evidence just by using transformations ( 10.13b 
and writing the corresponding spectrum in terms of the usual uj and k z spectral parameters. 

In the following, we consider some cases with /3 = and (3 > 0. 



0.2. 1. 1 Closed analytical expressions describing some ideal nondiffracting 
pulses 
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Let us first consider, in eq. tl0.15\) , the following spectra, of the type ( I6>.77l ) with f3 = 0: 

A(a,f3) = aV6{/3)e~ aVa (0.18) 

A(a,j3) = aVS(P)J (2dy/a)e- aVa (0.19) 

A(a,f3) = S(pf^^-e- aVa (0.20) 
a 

with a > and d being constants. 

One can obtain from the above spectra the following superluminal localized wave solu- 
tions, respectively: 

— From spectrum ( 10.181 1. we can use the identity (6.61 l.l)inref.[19], obtaining the well 
known ordinary X xave solution (also called X-shaped pulse) 

*(p,C) = X = aV (0.21) 

^(aV-*0 2 +(£-l) p 2 

— Using spectrum (10.19b and the identity (6.6444) of ref.[19], we get 



*0,C) - X-J [\j^--1 {aV)- 2 d 2 X 2 p ) cxp[-(aV -i()(aV)- 2 d 2 X 2 ] 



(0.22) 



— The superluminal nondiffracting pulse 

*(p,C) = sin -1 



4- r ( Jx-z + (d/aV) 2 + 2pd(aV)- 2 JV 2 /c 2 
aV \ v 

+ \jx- 2 + (d/aV) 2 - 2pd(aV)- 2 v / V 2 /c 2 \ 



(0.23) 



is obtained from spectrum (I0.201 > using identity (6.752.1) of ref.[19] for a > and d > 0. 

From the previous discussion, we get to know that any solution obtained from spectra 
of the type ( 10. 17b with (3q = is free from noncausal (backwards) components. 

In addition, when f3o = 0, we can see that the pulsed solutions depend on z and t through 
( = z — Vt only, and so propagate rigidly, i.e. without distortion. Such pulses can be 
transversally localized only if V > c, because if V = c the function ^ has to obey the 
Laplace equation on transverse planes[12,13]. 

Many others superluminal localized waves can be easily constructed[13] from the above 
solutions just by taking the derivatives (of any order) with respect to (. It is also possible 
to show[13] that the new solutions obtained in this way have their spectra shifted towards 
higher frequencies. 
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Now, let us consider, in eq. A0.15\) , a spectrum of the type ( I6>. 1 71 ) with /3q > 0: 

A(a,f3) = aVS(f3 - (3 a ) e - aVa (0.24) 

with a a positive constant. 

As we have seen, the presence of the delta function, with the constant (3q > 0, implies that 
we are integrating (summing) Bessel beams along the continuous line ui = Vk z + 2V/3q. 
Now, the function S(a) = aVexp(— aVuo) entails that we are considering a frequency 
spectrum of the type S(ui) oc cxp(— aoj), and therefore with a bandwidth given by Aw = 
1/a. 

Since /3o > 0, the interval Vfio < ui < 2Vfio (or, equivalently in this case, < a < f3o), 
corresponds to backward Bessel beams, i.e. negative values of k z . However, we can get 
physical solutions when making the contribution of this frequency interval negligible. In 
this case, it can be done by making aftoV << 1, so that the exponential decay of the 
spectrum S with respect to uo is very slow and the contribution of the interval u> > 2VPo 
(where k z > 0) overruns the V(3q < uj < 2Vfio (where k z < 0) contribution. 

Incidentally, we note that, once we ensure the causal behavior of the pulse by making 
aVPo << 1 in (10.24b . we have that Aa = 1/aV >> f3o, and so we can simplify the 
argument of the Bessel function, in the integrand of superposition (10.15b . by neglecting the 
term (V 2 /c 2 — 1)/3q- With this, the superposition (10.15b . with the spectrum (10.24b . can be 
written as 



(0.25) 

Now, we can use identity (6.6 16.1) of ref.[19] and obtain the new localized superluminal 
solution called[13] Superluminal Focus Wave Mode (SFWM): 



SFWM 



0, C, v) 



exp 



V 2 -c 2 



{{aV - iQ - aV X' 1 ) 



(0.26) 



where, as before, X is the ordinary X pulse ( 10.21b . The center of the SFWM is localized 
on p = and £ = (i.e. AT z = Vt). The intensity, l^l 2 , of this pulse propagates 
rigidly, it being a function of p and £ only. However, the complex function "Jsfwm (i.e. 
its real and imaginary parts) propagate just with local variations, recovering their whole 
three dimensional form after each space and time interval given by Azo = 7t//3q and 
Ai = tt/[3 V. 

The SFWM solution written above, for V — ► c + reduces to the well known Focus 
Wave Mode (FWM) solution [1 1], travelling with speed c: 



e -i0oV 

^fwm(p,C,v) = ac— — — exp 



ac — i( 



Pop 2 



ac — iQ 



(0.27) 
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Let us also emphasize that, since [3q > 0, the spectrum (10. 24b results to be constituted 
by angular frequenies ui > V(3q. Thus, our new solution can be used to construct high 
frequency pulses. 

0.2. 1.2 Finite energy nondiffracting pulses 

In this subsection, we will show how to get finite energy localized wave pulses. These 
new waves can propagate for long distances while maintaining their spatial resolution, i.e. 
they possess a large depth of field. 

As we have seen, ideal nondiffracting waves can be constructed by superposing Bessel 
beams Ceq. dO. 12b for cylindrical symmetry) with a spectrum A(oj, k z ) that satisfies a linear 
relationship between u> and k z . In the general bidirectional decomposition method, this can 
be obtained by using spectra of the type ( 10.171 ) in superposition ( 10.151 ). 

Solutions of this type possess an infinity depth of field, however they exhibit infinite 
energy[l 1,13]. To overcome this problem, we can truncate an ideal nondiffracting wave 
by a finite aperture, and the resulting pulse will have finite energy and a finite field depth. 
Even so, such field depths may be very large when compared with those of ordinary waves. 

The problem in this case is that the resulting field has to be calculated from the diffraction 
integrals (such as the well known Rayleigh-Sommerfeld formula) and, in general, a closed 
analytical formula for the resulting pulse cannot be obtained. 

However, there is another way to construct localized pulses with finite energy[13]. That 
is, by using spectra A(ui,k z ) in (10.12b whose domains are not restricted to be defined 
exactly over the straight line lu = Vk z + b, but around that line, where the spectra should 
concentrate their main values. In other words, the spectrum has to be well localized in the 
vicinity of that line. 

Similarly, in terms of the generalized bidirectional decomposition given in (10.15b . finite 
energy nondiffracting wave pulses can be constructed considering well localized spectral 
functions A(a, j3) in the vicinity of the line j3 = fto, being (3o a constant. 

To exemplify this method, let us consider the following spectrum 

( aqVe- aVa e- q( -' 3 - f3o '> for j3 > /3 

A(a,0) = \ (0.28) 
[ for < /3 < O 

in the superposition (10.15b . quantities a and q being free positive constants and V the peak's 
pulse velocity (here V > c). 

It is easy to see that the above spectrum is zero in the region above the (3 = (3q line, 
while it decays in the region below (as well as along) such a line. We can concentrate this 
spectrum on f3 = /3o by choosing values of q in such a way that qfto >> 1. The faster the 
spectrum decay takes place in the region below the (3 = (3q line, the larger the field depth 
of the corresponding pulse results to be. 
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Besides this, once we choose qj3$ >> 1 to obtain pulses with large field depth, we 
can also minimize the contribution of the noncausal (backward) components by choosing 
aVfto << 1; in analogy with the results we obtained for the SFWM case. 

Still in analogy with the SFWM case, when we choose qflo >> 1 (i.e. A long field 
depth) and aVfio < < 1 (minimal contribution of backward components), we can simplify 
the argument of the Bessel function, in the integrand of superposition ( 10.151 ). by neglecting 
the term (V 2 /c 2 - 1)$. 

With the observations above, we can write the superposition (10.15b with the spectrum 
( 10281 ) as 

*(px,v) « aqV S P d(3 1 daJo ( P y(5^0 a2 + 2 (^ + 1 ) 

X e-iPv^Cg-liP-M e ~ aVa 

(0.29) 

and, using identity (6.616.1) given in ref.[19], we get 



/>oo 



(0.30) 



which can be viewed as a superposition of the SFWM pulses (see eq. ( 10.261 )). 

The above integration can be easily made and results[13] in the so called Superluminal 
Modified Power Spectrum (SMPS) pulse: 

^to)=«I^M (0.31) 
q-{Y - vq) 

where X is the ordinary X pulse (10.211 ) and Y is defined as 

Y = v* 2 -* ^ aV ~ zC) ~ aVX ^ (032) 

The SMPS pulse is a superluminal localized wave, with field concentration around p = 
and C = (i.e. in z — Vt), and with finite total energy. We will show that the depth of 
field, Z, of this pulse is given by Zsmps = ?/2. 

An interesting property of the SMPS pulse is related with its transverse width (the 
transverse spot size at the pulse center). It can be shown from (10.31b that for the cases 
where aV << 1/(3q and q/3a >> 1, i.e., for the cases considered by us, the transverse spot 
size, Ap, of the pulse center (£ = 0) is dictated by the exponential function in (10.31b and is 
given by 



aV V 2 - c 2 

Af> = C I/ p (V 2 + c 2 ) + 4[3 2 (V 2 + c 2 ) 2 (033) 



X 



which clearly does not depend on z, and so remains constant during the propagation. In 
other words, in spite of the fact that the SMPS pulse suffers an intensity decrease during 
the propagation, it preserves its transverse spot size. This interesting characteristic is not 
verified in ordinary pulses, like the gaussian ones, where the amplitude of the pulse decreases 
and the width increases by the same factor. 

Figure 10. l.l shows a SMPS pulse intensity, with (3 = 33 m -1 , V = 1.01c, a = 10~ 12 s 
and q — 10 5 m, at two different moments, for t — and after 50 km of propagation, where, 
as we can see, the pulse becomes less intense (half of its initial peak intensity). It can be 
noted that, in spite of the intensity decrease, the pulse maintains its transverse width, as one 
can see from the 2D plots in Fig. ( 10.1. 1 ), which show the field intensities in the transverse 
sections at z = and z = q/2 = 50 km. 





Figure 0.1. Representation of a Superluminal Modified Power Spectrum pulse, eq. dO.3 It . Its total 
energy is finite (even without any truncation), and so it gets deformed while propagating, since its 
amplitude decreases with time. In Fig lO. l.h we represent, for t = 0, the pulse corresponding to 
Po — 33 m -1 , V — 1.01c, a — 10~ 12 s and q — 10 5 m. In Fig lO. l.b it is depicted the same pulse 
after having travelled 50 km. 

Other three important well known finite energy nondiffracting solutions can be obtained 
directly from the SMPS pulse: 

— The first one, obtained from ( 10.311 1 by making /3q = 0, is the so called[13] Superlu- 
minal Splash Pulse (SSP), 

*ssp(p,C,'/) = J X v (034) 
q + if] — Y 

— The other two are luminal pulses. By taking the limit V — > c + in the SMPS pulse 
d0.3U . we get the well known[l 1] luminal Modified Power Spectrum (MPS) pulse 

*mps(p,C,?7 = / , ■ w ..x , 2 exp (0.35) 

[q + iT])(ac — iQ + p z \ac — iQJ 

Finally, by taking the limit V — > c + and making (3q — in the SMPS pulse (or, 
equivalently, by making (3 = in the MPS pulse ( 10.35b . or,instead, by taking the limit 
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V -> c+ in the SSP (10.34b ). we obtain the well knownfll] luminal Splash Pulse (SP) 
solution 

*sp(p,C,»7) = 7 - ■ w° gC 2 (°- 36 ) 

(q + ir])(ac - iQ + p z 

It is also interesting to notice that the X and SFWM pulses can be obtained from the SSP 

and SMPS pulses (respectively) by making q — > oo in Eqs. (10.34l) and (10.311) . As a matter 

of fact, the solutions SSP and SMPS can be viewed as THE finite energy versions of the X 

and SFWM pulses, respectively. 



Some characteristics of the SMPS pulse: 



Let us examine the on-axis (p = 0) behavior of the SMPS pulse. 
On p = we have 

*smpsG° = 0,C,?7) = aqVe-^KaV -iOiq + ivT 1 (0.37) 

From this expression, we can show that the longitudinal localization Az, for t — 0, of 
the SMPS pulse square magnitude is 



Az = 2aV (0.38) 

If we now define the field depth Z as the distance over which the pulse's peak intensity 
is 50% at least of its initial valueQ, then we can obtain, from ( 10.37b . the depth of field 

Zsmps = | (0.39) 

which depends only on q, as we expected since q regulates the concentration of the spectrum 
around the line u> — Vk z + 2Vpo- 

Now, let us examine the maximum amplitude M of the real part of (10.37b , which for 
z = Vt writes = and r\ = 2z): 

^ r T / . cos(28 z) — 2(z/q)sm(28 z) ,„ 

M SM ps = Re*sMPs(p = 0,z = V r f) = , 77 ,L (0-40) 

1 + ^{z/qY 

Initially, for z = 0, t = 0, one has M = 1 and can also infer that: 
(i) when z/q « 1, namely, when z << Z, equation (10.40b becomes 

M SM ps ~ cos(2/3 2;) for z « Z (0.41) 

and the pulse's peak actually oscillates harmonically with "wavelength" Azo = tt/Bq and 
"period" A^o = k/VPo, all along its field depth. 

* We can expect that while the pulse peak intensity is maintained, also is its spatial form. 
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(ii) When z/q » 1, namely z >> Z, equation ( 10.401 ) becomes 

M SMPS «-^M f 0rz »Z (0.42) 
2z/q 

Therefore, beyond its depth of field, the pulse goes on oscillating with the same Azo, 
but its maximum amplitude decays proportionally to z. 

In the next two Sections we are going to see an interesting application of the localized 
wave pulses. 



0.3 SPACE-TIME FOCUSING OF X-SHAPED PULSES 



In this Section we are going to show how one can in general use any known Superluminal 
solution, to obtain from it a large number of analytic expressions for space-time focused 
waves, endowed with a very strong intensity peak at the desired location. 

The method presented here is a natural extension of that developed by A. Shaarawi et 
al.[20], where the space-time focusing was achieved by superimposing a discrete number 
of ordinary X-waves, characterized by different values 9 of the axicon angle. 

In this section, based on ref.[21], we will go on to more efficient superpositions for 
varying velocities V, related to 8 through the known[3,4] relation V = cj cosO. This 
enhanced focusing scheme has the advantage of yielding analytic (closed-form) expressions 
for the spatio-temporally focused pulses. 

Let us start considering an axially symmetric ideal nondiffracting superluminal wave 
pulse ip(p,z — Vt) in a dispersionless medium, where V = cj cos 6 > c is the pulse 
velocity, 9 being the axicon angle. As we have seen in the previous Section, pulses like 
these can be obtained by a suitable frequency superposition of Bessel beams. 

Suppose that we have now N waves of the type ip n (p, z — Vn(t — tn)), with different 
velocities c < V\ < V2 < .. < Vn, and emitted at (different) times t n ; quantities t n being 
constants, while n = 1, 2, ...N. The center of each pulse is located at z = V n (t — t n ). To 
obtain a highly focused wave, we need all the wave components ipn(p, z — Vn(t — tn)) to 
reach the given point, z = Z{, at the same time t = t[. On choosing t\ = for the slowest 
pulse ipi, it is easily seen that the peak of this pulse reaches the point z = Zf at the time 
tf = Z{/Vi. So we obtain that, for each ip n , the instant of emission t n must be 

(0.43, 

With this, we can construct other exact solutions to the wave equation, given by 



*(p, *,t) = J v A(V) i> (p, z - V (t - (-^- - 7) * ) ) » (0-44) 

where V is the velocity of the wave ip(p, z—Vt) in the integrand of (10.44t . In the integration, 
V is considered as a continuous variable in the interval [V m i n , V max ]. In eq. (2), A(V) is 
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the velocity-distribution function that specifies the contribution of each wave component 
(with velocity V) to the integration. The resulting wave $f(p, z, t) can have a more or 
less strong amplitude peak at z — Z{, at time if = Zf/V m i n , depending on A(V) and on 
the difference V max — V m i n . Let us notice that also the resulting wavefield will propagate 
with a Superluminalpeak velocity, depending on A(V) too. In the cases when the velocity- 
distribution function is well concentrated around a certain velocity value, one can expect the 
wave ( I0.441 > to increase its magnitude and spatial localization while propagating. Finally, 
the pulse peak acquires its maximum amplitude and localization at the chosen point z = Z{, 
and at time t = Z{ / V m i n , as we know. Afterwards, the wave suffers a progressive spreading, 
and a decreasing of its amplitude. 



0.3.1 Focusing Effects by Using Ordinary X-Waves 

Here, we present a specific example by integrating (I0.44t over the standard, ordinary [4] 
X-waves,X = aV[(aV- i{z - Vt)) 2 + (V 2 /c 2 - Vjp 2 ]- 1 ' 2 . When using this ordinary X- 
wave, the largest spectral amplitudes are obtained for low frequencies. For this reason, one 
may expect that the solutions considered below will be suitable mainly for low frequency 
applications. 

Let us choose, then, the function tjj in the integrand of eq. (l0.44b to be ifj(p,z,t) = 
X(p, z-V(t- (1/tU. - 1/V)zt)), viz. 



i>(p,z,t) = x = 



iV 



After some manipulations, one obtains the analytic integral solution 



with 



V(p,z,t) 



P 



v maJ 



aVA(V) 



lin ^JPV 2 +QV + R 
(t - jr- ~ ai) {z ~ Zf) 



-AV 



Q = 2 

R = [-{z-z t ) 2 -p 2 



(0.45) 



(0.46) 



(0.47) 



In what follows, we illustrate the behavior of some new spatio-temporally focused pulses, 
by taking into consideration some different velocity distributions A(V). These new pulses 
are closed analytical exact solutions of the wave equation. 



First example: 
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Let us consider our integral solution ( 10.46t with A(V) = ls/m. In this case, the 
contribution of the X-waves is the same for all velocities in the allowed range [Vmin, VWx]- 
Using identity 2.264.2 listed in ref.[19], we get the particular solution 




(0.48) 



where P, Q and R are given in eq.( l0.47t . A 3-dimensional (3D) plot of this function is 
provided in Fig |D~2~l where we have chosen a = 10~ 12 s, V m m = 1-001 c, V ma ^ = 1.005 c 
and Z[ = 200 cm. It can be seen that this solution exhibits a rather evident space-time 
focusing. An initially spread-out pulse (shown for t — 0) becomes highly localized at 
t = U — Z[/V m in = 6.66 ns, the pulse peak amplitude at Z{ being 40.82 times greater than 
the initial one. In addition, at the focusing time tf the field is much more localized than at 
any other times. The velocity of this pulse is approximately V = 1.003 c. 




Figure 0.2. Space-time evolution of the Superluminal pulse represented by eq. d0.48l l: the chosen 
parameter values are a — 10~ 12 s; Vxmn = 1.001 c; Knax = 1.005 c while the focusing point 
is at Zf = 200 cm. One can see that this solution is associated with a rather good spatio-temporal 
focusing. The field amplitude at z = zt is 40.82 times larger than the initial one. The field amplitude 
is normalized at the space-time point p = 0, z = Zf , t — tf. 



Second example: 



In this case we choose A(V) — 1/V (s/m), and, using the identity 2.261 in ref.[19], 
eq. ( 10.46b gives 
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n t) = a ^ ( 2V P(PV^ + QV^ + -R) + 2PV^ + Q\ 
VP \ 2 y/P^PV^ + QV min + R) + 2PV min + Q J 

Other exact closed-form solutions can be obtained[2 1 ] considering, for instance, velocity 
distributions like A(V) = l/V 2 and A(V) = 1/V 3 . 

Actually, we can construct many others spatio-temporally focused pulses from the above 
solutions, just by taking their time derivatives (of any order). It is also possible to show[21] 
that the new solutions obtained in this way have their spectra shifted towards higher fre- 
quencies. 



0.4 CHIRPED OPTICAL X-TYPE PULSES IN MATERIAL MEDIA 

The theory of the localized waves was initially developed for free space (vacuum). In 
1996, Sonajalg et al.[5] showed that the localized wave theory can be extended to include 
(unbounded) dispersive media. This was obtained by making the axicon angle of the 
Bessel beams (BBs) vary with the frequency[5-7] in such a way that a suitable frequency 
superposition of these beams does compensate for the material dispersion. Soon after 
this idea was reported, many interesting nondiffracting/nondispersive pulses were obtained 
theoretically[5-7] and experimentally [5]. 

In spite of this extended method to be of remarkable importance, working well in theory, 
its experimental implementation is not so simpl J+J. 

In 2004 Zamboni-Rached et al.[22] developed a simpler way to obtain pulses capable of 
recovering their spatial shape, both transversally and longitudinally, after some propagation. 
It consisted in using chirped optical X-typed pulses, while keeping the axicon angle fixed. 
Let us recall that, by contrast, chirped Gaussian pulses in unbounded material media may 
recover only their longitudinal shape, since they undergo a progressive transverse spreading 
while propagating. 

The present section is devoted to this approach. 

Let us start with an axis-symmetric Bessel beam in a material medium with refractive 
index n(u>): 

tp(p,Z,t) = Jo(fcpp) exp(i/3z) exp(— iu)t) , (0.50) 

where it must be obeyed the condition k 2 = n 2 (lu)ui 2 / c 2 — (3 2 , which connects among 
themselves the transverse and longitudinal wave numbers k p and /3, and the angular fre- 
quency uj. In addition, we impose that k 2 > and uj/(3 > 0, to avoid a nonphysical 
behavior of the Bessel function Jq(.) and to confine ourselves to forward propagation only. 

tWe refer the interested reader to quotations [5-7] for obtaining a description, theoretical and experimental, of 
this extended method 
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Once the conditions above are satisfied, we have the liberty of writing the longitudinal 
wave number as (3 — (n(u>)uj cos 9) jc and, therefore, k p = (n(w)ui sin 9) / c; where (as in 
the free space case) 8 is the axicon angle of the Bessel beam. 

Now we can obtain a X-shaped pulse by performing a frequency superposition of these 
Bessel beams [BB], with /3 and k p given by the previous relations: 



S(lu) Jo ( — sintfpj exp[iP(u)z] exp(— iujt) duj , 



(0.51) 



where S(u>) is the frequency spectrum, and the axicon angle is kept constant. 

One can see that the phase velocity of each BB in our superposition ( 10.5 U is different, 
and given by V^hasc = c/(n(w) cos 9). So, the pulse represented by eq. (10. 5 11 1 will suffer 
dispersion during its propagation. 

As we said, the method developed by Sonajalg et al.[5] and explored by others[6,7], 
to overcome this problem, consisted in regarding the axicon angle 9 as a function of the 
frequency, in order to obtain a linear relationship between (3 and to. 

Here, however, we wish to work with a fixed axicon angle, and we have to find out 
another way for avoiding dispersion and diffraction along a certain propagation distance. 
To do that, we might choose a chirped gaussian spectrum S(ui) in eq. (l0.51b 

S(u) = fo exphgV - c ) 2 ] with q 2 = ^ , (0.52) 

V2tt(1 + iC) 2(1 + iG) 

where is the central frequency of the spectrum, To is a constant related with the initial 
temporal width, and C is the chirp parameter (we chose as temporal width the half-width of 
the relevant gaussian curve when its heigth equals 1/e times its full heigth). Unfortunately, 
there is no analytical solution to eq. (10. 5 11 1 with S(uS) given by eq. (10.52b . so that some 
approximations are to be made. 

Then, let us assume that the spectrum S{uS), in the surrounding of the carrier frequency 
u>o , is enough narrow that Auj/ujq << 1, so to ensure that (3{lo) can be approximated by 
the first three terms of its Taylor expansion in the vicinity of ojq: That is, (3(uj) « + 
/3'(w)\u (w — ujo) + (1/2)0" (u))\ Uo (uj — wo) 2 ; where, after using (3 — n(uj)uj cos 6>/c, it 
results that 



d(3 cos 9 
8lo c 



dn 



d 2 fi cos( 
dui 2 c 



^ dn d 2 n 
du> duj 2 



(0.53) 



As we know, (3 1 (u>) is related to the pulse group-velocity by the relation Vg — l//3'(u>). 
Here we can see the difference between the group-velocity of the X-type pulse (with a fixed 
axicon angle) and that of a standard gaussian pulse. Such a difference is due to the factor 
cos 9 in eq. ( 10.53b . Because of it, the group-velocity of our X-type pulse is always greater 
than the gaussian's. In other words, (V g )x = (1/ cos 9) (V g ) ga uss- 
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We also know that the second derivative of(3(u>) is related to the group-velocity dispersion 
(GVD) /3 2 by /3 2 = (3"{lo). 

The GVD is responsible for the temporal (longitudinal) spreading of the pulse. Here 
one can see that the GVD of the X-type pulse is always smaller than that of the standard 
gaussian pulses, due the factor cos 9 in eq. (10.53b . Namely: (/3 2 )x = cos 0(j3 2 ) gajVLSS . 

Using the above results, we can write 



V(p,z,t) = 



To exp[if3(u>o) z ] exp(— iwot) 

y/2TT(l+iC) 



'n(u>)w . 
dui Jq ( sin p 



x exp ji — [z — V g t]j exp |(w — uoY 



102 2 
— z — a 
2 q 



(0.54) 



The integral in eq. ( 10.54b cannot be solved analytically, but it is enough for us to obtain 
the pulse behavior. Let us analyze the pulse for p — 0. In this case we obtain: 



T exp[if3(u)o)z] exp(— iu>ot) 
^T 2 -i(3 2 (l + iCYz 



exp 



-(z-V E i) 2 {l + iC) 



2V?[T 2 -H3 2 {l + iC)z] 



(0.55) 

From eq. (10.55b one can immediately see that the initial temporal width of the pulse 
intensity is Tq and that, after a propagation distance z, the time-width T% becomes 



T 



1 



C(3 2 z 
T 2 



T 2 



1/2 



(0.56) 



Relation ( 10.56b describes the pulse spreading-behavior. One can easily show that such 
a behavior depends on the sign (positive or negative) of the product fcp, as is well known 
for the standard gaussian pulses[23]. 

In the case f3 2 C > 0, the pulse will monotonically become broader and broader with 
the distance z. On the other hand, if (3 2 C < the pulse will suffer, in a first stage, a 
narrowing, and then it will spread during the rest of its propagation. So, there will be 
a certain propagation distance AT which the pulse will recover its initial temporal width 
(7\ = T ). From relation ( 10.561 ). we can find this distance Zti=t (considering [3 2 C < 0) 
to be 



Zt-i =Tr, — 



-2CT 2 



(0.57) 



fc(C 2 + l) ' 

One may notice that the maximum distance at which our chirped pulse, with given To 
and f3 2 , may recover its initial temporal width can be easily evaluated from eq. (10.57b . and 
it results to be Ldi sp = Tq /(3 2 . We shall call such a maximum value Ldisp the "dispersion 
length". It is the maximum distance the X-type pulse may travel while recovering its initial 
longitudinal shape. Obviously, if we want the pulse to reassume its longitudinal shape at 
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some desired distance z < Tdisp, we have just to suitably choose the value of the chirp 
parameter. 

Let us emphasize that the property of recovering its own initial temporal (or longitu- 
dinal) width may be verified to exist also in the case of chirped standard gaussian pulses. 
However, the latter will suffer a progressive transverse spreading, which will not be re- 
versible. The distance at which a gaussian pulse doubles its initial transverse width wo 
is Zdiff = V3kwq/ Ao, where Ao is the carrier wavelength. Thus, we can see that optical 
gaussian pulses with great transverse localization will get spoiled in a few centimeters or 
even less. 

Now we shall show that it is possible to recover also the transverse shape of the chirped 
X-type pulse intensity; actually, it is possible to recover its entire spatial shape after a 
distance Zt 1= t - 

To see this, let us go back to our integral solution ( 10.541 i. and perform the change of 
coordinates (z, t) — ► (Az, t c = z c /V g ), with 



z = z c + Az 



t = t c = — 



(0.58) 



where z c is the center of the pulse (Az is the distance from such a point), and t c is the 
time at which the pulse center is located at z c . What we are going to do is comparing 
our integral solution (10.541 i. when z c = (initial pulse), with that when z c = Zt 1= t = 
-2CT*/(f3 2 (C 2 + 1)). 

In this way, the solution (10.54b can be written, when z c = 0, as 



t / „ . . T exp(i/3 Az) , , 

®{p,Zc = 0,Az) = " == - / &>J Q (k p (u>)p) exp 



-T 2 (to - luo) 2 



2(1 + C 2 ) 



x exp < i 



(o;-w )Az (lo-loq) 2 /3 2 Az (w-w )T 2 C 
V e + 2 + 2(1 + C 2 ) 



(0.59) 

where we have taken the value q given by ( 10.52l >. 

To verify that the pulse intensity recovers its entire original form at z c = Zt 1= t = 
—2 CTq l\fii(C 2 + 1)], we can analyze our integral solution at that point, obtaining 
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T exp <^ ifio 



z c = Z Ti =t ,Az) 



-Az'- 



cos0n(cJo)V g 



^(l + iC) 



x / dw Jo(k p ((jj)p) exp 



-T 2 (u - uj ) 2 



2(1 + C 2 ) 



(0.60) 



x exp < —i 



(oj-m a )Az' (uj-u Ja ) 2 p 2 Az / {lo-lo ) 2 T 2 C 
V g 2 2(1 + C 2 ) 



where we put Az = —Az'. In this way, one immediately sees that 



z c = 0,Az) 



|*09, z c = Z Tl=To ,-Az) 



(0.61) 



Therefore, from eq. (l0.6U it is clear that the chirped optical X-type pulse intensity will 
recover its original three-dimensional form, with just a longitudinal inversion at the pulse 
center: The present method being, in this way, a simple and effective procedure for com- 
pensating the effects of diffraction and dispersion in an unbounded material medium; and 
a method simpler than the one of varying the axicon angle with the frequency. 

Let us stress that we can choose the distance z = Zt 1= t < ^disp at which the pulse 
will take on again its spatial shape by choosing a suitable value of the chirp parameter. 

Till now, we have shown that the chirped X-type pulse recovers its three-dimensional 
shape after some distance, and we have also obtained an analytic description of the pulse 
longitudinal behavior (for p = 0) during propagation, by means of eq. (10.55b . However, 
one does not get the same information about the pulse transverse behavior: We just know 
that it is recovered at z = Zt 1= t - 

So, to complete the picture, it would be interesting if we could get also the transverse 
behavior in the plane of the pulse center z — V g t. In that way, we would obtain quantitative 
information about the evolution of the pulse-shape during its entire propagation. 

We will not examine the mathematical details here; we just affirm that THE transverse 
behavior of the pulse (in the plane z = z c = V g t) during its whole propagation can 
approximately be described by 



XX 



V(p,z= z c ,t= z c /V g ) 



T exp[i/3(wo)z] exp(— iuj^t) 

y/2-K{l+iC) 



exp 



tan 2 p 2 



8V 2 (-i(3 2 z c /2 + q 2 )_ 



+ q 2 



r(i/2)j o r M " sin ^ ')/o 



tan 2 p 2 



8V?(-il3 2 z c /2 +q 2 ) 



+2 



E 

P =i 



2PT( P +l/2)T(p+l) 

r(2 P + i) 



2 1, 



n(u>o)ujo sin 9 p 



»2p 



tan 2 p 2 



8K 2 (-iftz c /2 + g 2 ) 



(0.62) 



where I p (.) is the modified Bessel function of the first kind of order p, quantity T(.) being 
the gamma function and q given by ( 10.52b . 

The interested reader can consult ref.[22] for details on how to obtain eq. ( 10.62b from 

At a first sight, this solution could appear very complicated, but the series in its r.h.s. 
gives a negligible contribution. This fact renders our solution ( 10.62b of important practical 
interest and we will use it in the following. 

For additional information about the transverse pulse evolution (extracted from eq. (10.62b ). 
the reader can consult again ref.[22] . In the same reference, it is analyzed the effect of a 
finite aperture generation on the chirped X-type pulses. 



0.4.1 An example: Chirped Optical X-typed pulse in bulk fused Silica 



For a bulk fused Silica, the refractive index n(w) can be approximated by the Sellmeier 
equation[23] 



N B J2 

n 2 (cj) = 1 + V 2 J \ , (0.63) 

3=1 3 

where ujj are the resonance frequencies, Bj the strength of the jth resonance, and N the 
total number of the material resonances that appear in the frequency range of interest. For 
our purposes it is appropriate to choose N = 3, which yields, for bulk fused silica[23], the 
values B 1 = 0.6961663; B 2 = 0.4079426; B 3 = 0.8974794; Ai = 0.0684043 p m; 
A 2 = 0.1162414 p m; and A 3 = 9.896161 p m. 
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Now, let us consider in this medium a chirped X-type pulse, with Ao = 0.2 /im, To = 
0.4 ps, C = — 1 and with an axicon angle 6 = 0.00084 rad, which correspond to an initial 
central spot with Apo = 0.117 mm. 

We get from eqs. (10.55b and (10.621 1. the longitudinal and transverse pulse evolution, 
which are represented by Fig l0.3.l 




Figure 0.3. (a): Longitudinal-shape evolution of a chirped X-type pulse, propagating in fused silica 
with Ao = 0.2// m, To = 0.4 ps, C = — 1 and axicon angle 8 = 0.00084 rad, which correspond 
to an initial transverse width of Apo = 0.117 mm. (b): Transverse-shape evolution for the same 
pulse. 

From Fig |Q.3.f a). we can notice that, initially, the pulse suffers a longitudinal narrowing 
with an increase of intensity till the position z = Tq /2/?2 = 0.186m. After this point, the 
pulse starts to broaden decreasing its intensity and recovering its entire longitudinal shape 
(width and intensity) at the point z = Tq 2 //?2 = 0.373m, as it was predicted. 

At the same time, from Fig |Q.3.f b), one can notice that the pulse maintains its transverse 
width Ap = 2.4 c/(n(wo)wo s'm9) = 0.117 mm (because Tqljq >> 1) during its entire 
propagation; however, the same does not occur with the pulse intensity. Initially the pulse 
suffers an increase of intensity till the position z c — Tq/2/3 2 = 0.186m; after this point, 
the intensity starts to decrease, and the pulse recovers its entire transverse shape at the point 
z c = Tq / P2 = 0.373m, as it was expected by us. Here we have skipped the series on the 
r.h.s. of eq. (10.621 i. because, as we already said, it yields a negligible contribution. 

Summarizing, from Fig l0.3.l we can see that the chirped X-type pulse recovers totally 
its longitudinal and transverse shape at the position z — Tnisp — Tq / (3 2 = 0.373 m, as we 
expected. 

Let us recall that a chirped gaussian pulse may just recover its longitudinal width, but 
with an intensity decrease, at the position given by z — Zt 1= t — ^disp = Tq//32. Its 
transverse width, on the other hand, suffers a progressive and irreversible spreading. 
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0.5 MODELING THE SHAPE OF STATIONARY WAVE FIELDS: FROZEN 
WAVES 

In this Section we develop a very simple method[10,16,17], by having recourse to su- 
perpositions of forward propagating and equal-frequency Bessel beams, that allows one 
controlling the longitudinal beam-intensity shape within a chosen interval < z < L, 
where z is the propagation axis and L can be much greater than the wavelength A of the 
monochromatic light (or sound) which is being used. Inside such a space interval, indeed, 
we succeed in constructing a stationary envelope whose longitudinal intensity pattern can 
approximately assume any desired shape, including, for instance, one or more high-intensity 
peaks (with distances between them much larger than A); and which results — in addition — 
to be naturally endowed also with a good transverse localization. Since the intensity en- 
velopes remains static, i.e. with velocity V = 0, we call "Frozen Waves" (FW) such new 
solutions [10, 16, 17] to the wave equations. 

Although we are dealing here with exact solutions of the scalar wave equation, vectorial 
solutions of the same kind for the electromagnetic field can be obtained, since solutions to 
Maxwell's equations follow naturally from the scalar wave equation solutions[24,25]. 

First, we present the method considering lossless media[16,17] and, in the second part 
of this section, we extend the method to absorbing media[10]. 

0.5.1 Stationary wavefields with arbitrary longitudinal shape in lossless 
media, obtained by superposing equal-frequency Bessel beams 

Let us start from the well-known axis-symmetric zeroth order Bessel beam solution to the 
wave equation: 

1>(p,z,t) = MkpPieW'e-*"* (0.64) 

with 

k 2 p = ^ - 1 , (0.65) 

where u, k p and (3 are the angular frequency, the transverse and the longitudinal wave 
numbers, respectively. We also impose the conditions 

lu/P > and k 2 p > (0.66) 

(which imply u/(3 > c) to ensure forward propagation only (with no evanescent waves), as 
well as a physical behavior of the Bessel function J . 

Now, let us make a superposition of 2N + 1 Bessel beams with the same frequency uj , 
but with different (and still unknown) longitudinal wave numbers fi m : 

N 

*(p,z,t) = e~ iuot A ™MKmP)e il3 ™ z , (0.67) 

m=-N 
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where thye m represent integer numbers and the A m are constant coefficients. For each 
m, the parameters ui , k pm and (3 m must satisfy ( I0.65l l. and, because of conditions ( 10.66b . 
when considering luq > 0, we must have 

< p m < — ■ (0.68) 

c 

Let us now suppose that we wish \^(p, z, t)\ 2 , given by eq.( l0.67l i. to assume on the axis 
p = the pattern represented by a function |_F(z)| 2 , inside the chosen interval < z < L. 
In this case, the function F(z) can be expanded, as usual, in a Fourier series: 



F(z) = J2 B ™ e ^ 

m— — oo 

where 



B m = ijf F(z)e- l ^ mz dz . 

More precisely, our goal is finding out, now, the values of the longitudinal wave numbers 
P m and the coefficients A m of (10.67b . in order to reproduce approximately, within the said 
interval < z < L (for p = 0), the predetermined longitudinal intensity-pattern \F(z) | 2 . 
Namely, we wish to have 

2 

w \F(z)\ 2 with < z < L . (0.69) 

Looking at eq. ( 10.691 1. one might be tempted to take f3 m = 2irm/L, thus obtaining a 
truncated Fourier series, expected to represent approximately the desired pattern F(z). 
Superpositions of Bessel beams with f3 m = 2nm / L have been actually used in some works 
to obtain a large set of transverse amplitude profiles [26]. However, for our purposes, this 
choice is not appropriate, due to two principal reasons: 1) It yields negative values for 
f3 m (when m < 0), which implies backwards propagating components (since ujq > 0); 2) 
In the cases when L >> Ao, which are of our interest here, the main terms of the series 
would correspond to very small values of /3 m , which results in a very short field-depth of the 
corresponding Bessel beams (when generated by finite apertures), preventing the creation 
of the desired envelopes far form the source. 

Therefore, we need to make a better choice for the values of (3 m , which permits forward 
propagation components only, and a good depth of field. This problem can be solved by 
putting 

0m = Q + ^-m, (0.70) 



N 
m=-N 
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where Q > is a value to be chosen (as we shall see) according to the given experimental 
situation, and the desired degree of transverse field localization. Due to eq. ( 10.681 ). we get 

„ 2 7T LUn 

0<Q± — N< — . (0.71) 

L c 

Inequality ( 10.7 11 1, can be used to determine the maximum value of m, that we call N mSLX , 
once Q, L and ujq have been chosen. 

As a consequence, for getting a longitudinal intensity pattern approximately equal to the 
desired one, |F(z)| 2 , in the interval < z < L, eq. (10.67b should be rewritten as 

N 

V(p = 0,z,t) = e - lUot e lQz A m e z2 ? mz , (0.72) 

m=—N 

with 

A m = I / F(z)e- iS P m *dz. (0.73) 
L Jo 

Obviously, one obtains only an approximation to the desired longitudinal pattern, because 
the trigonometric series ( 10.72b is necessarily truncated (N < N max ). Its total number of 
terms, let us repeat, will be fixed once the values of Q, L and uq are chosen. 

When p 7^ 0, the wavefield fy(p,z, t) becomes 

N 

9{p,z,t) = e- iUot e^ z A m J (k pmP )e^ mz , (0.74) 

m=—N 

with 

k 2 pm = c 2 -(q + ^) 2 . (0.75) 

The coefficients A m will yield the amplitudes and the relative phases of each Bessel 
beam in the superposition. 

Because we are adding together zero-order Bessel functions, we can expect a high field 
concentration around p = 0. Moreover, due to the known non-diffractive behavior of the 
Bessel beams, we expect that the resulting wavefield will preserve its transverse pattern in 
the entire interval < z < L. 

The methodology developed here deals with the longitudinal intensity pattern control. 
Obviously, we cannot get a total 3D control, due the fact that the field must obey the wave 
equation. However, we can use two ways to have some control over the transverse behavior 
too. The first is through the parameter Q of eq. (10.70b . Actually, we have some freedom in 
the choice of this parameter, and FWs representing the same longitudinal intensity pattern 
can possess different values of Q. The important point is that, in superposition (10.74b . using 
a smaller value of Q makes the Bessel beams possess a higher transverse concentration (be- 
cause, on decreasing the value of Q, one increases the value of the Bessel beams transverse 
wave numbers), and this will reflect in the resulting field, which will present a narrower 
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central transverse spot. The second way to control the transverse intensity pattern is using 
higher order Bessel beams, and we shall show this in Section 5.1.1. 
Now, let us present a few examples of our methodology. 



First example: 



Let us suppose that we wish an optical wavefield with Ao = 0.632 /im, i.e. with 
ojo = 2.98 x 10 15 Hz, whose longitudinal pattern (along its z-axis) in the range < z < L 
is given by the function 



F(z) 



(z-h)(z-l 2 ) 
(h-h) 2 

1 

(z - l 5 )(z - l 6 ) 

{k-h? 





for h < z < l 2 
for I3 < z < Z4 
for h < z < Iq 



(0.76) 



elsewhere , 

where Z x = L/5 — Az 12 and l 2 = L/5 + Az 12 with Az 12 = L/50; while l 3 = 
L/2- Az 34 and l 4 = L/2 + Az 34 with Az 34 = L/10; and, at last, l 5 = 4L/5 - Az 56 
and 1% = 4L/5 + AZ56 with AZ56 = L/50. In other words, the desired longitudinal shape, 
in the range < z < L, is a parabolic function for h < z < l 2 , a unitary step function for 
I3 < z < I4, and again a parabola in the interval I5 < z < 1$, being zero elsewhere (within 
the interval < z < L, as we said). In this example, let us put L = 0.2 m. 

We can then easily calculate the coefficients A m , which appear in the superposition 
( 10.741 ). by inserting eq.( 10.76t into eq. ( 10.731 1. Let us choose, for instance, Q — 0.999 loq/c. 
This choice permits the maximum value N niax = 316form, as one can infer from eq. d0. 7 U . 
Let us emphasize that one is not compelled to use just TV = 316, but can adopt for N any 
values smaller than it; more generally, any value smaller than that calculated via inequality 
( 10.711 ). Of course, WHEN using the maximum value allowed for N, one gets a better result. 

In the present case, let us adopt the value N = 30. In Fig lOAf a) we compare the 
intensity of the desired longitudinal function F(z) with that of the Frozen Wave, ^(p = 
0, z, t), obtained from eq. (10.721) by adopting the mentioned value TV = 30. 



One can verify that a good agreement between the desired longitudinal behavior and 
our approximate FW is already got with TV = 30. The use of higher values for TV can only 
improve the approximation. Figure 10.4. l b) shows the 3D-intensity of our FW, given by 
eq. (10.741 ). One can observe that this field possesses the desired longitudinal pattern, while 
being endowed with a good transverse localization. 



Second example (controlling the transverse shape too): 
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m 2 



m 2 




p(mm) 



z(m) 



z(m) 



(a) 



(b) 



Figure 0.4. (a) Comparison between the intensity of the desired longitudinal function F(z) and 
that of our Frozen Wave (FW), Sb(p = 0, z, t), obtained from eq. jQ.72b . The solid line represents 
the function F(z), and the dotted one our FW. (b) 3D-plot of the field-intensity of the FW chosen in 
this case by us. 

We wish to take advantage of this example for addressing an important question: We 
can expect that, for a desired longitudinal pattern of the field intensity, by choosing smaller 
values of the parameter Q one will get FWs with narrower transverse width [for the same 
number of terms in the series entering eq. (l0.74H . because of the fact that the Bessel beams in 
eq. (I0.741 i will possess larger transverse wave numbers, and, consequently, higher transverse 
concentrations. We can verify this expectation by considering, for instance, inside the usual 
range < z < L, the longitudinal pattern represented by the function 



with l\ = L/2 — Az and I2 — L/2 + Az. Such a function has a parabolic shape, with its 
peak centered at L/2 and with longitudinal width 2Az/^/2. By adopting Ao = 0.632 /im 
(that is, luq = 2.98 x 10 15 Hz), let us use the superposition ( 10.74t with two different values 
of Q: We shall obtain two different FWs that, in spite of having the same longitudinal 
intensity pattern, will possess different transverse localizations. Namely, let us consider 
L = 0.06 m and Az = L/100, and the two values Q = 0.999 lo /c and Q = 0.995 lu /c. 
In both cases the coefficients A m will be the same, calculated from eq. (10.731 >. using this 
time the value N = 45 in superposition (10.74b . The results are shown in Figs. 3(a) and 3(b). 
Both FWs have the same longitudinal intensity pattern, but the one with the smaller Q is 
endowed with a narrower transverse width. 

With this, we can get some control on the transverse spot size through the parameter Q. 
Actually, eq. (10.74b . which defines our FW, is a superposition of zero-order Bessel beams, 
and, due to this fact, the resulting field is expected to possess a transverse localization 
around p = 0. Each Bessel beam in superposition ( 10.741 ) is associated with a central spot 




for h < z < h 



elsewhere 



(0.77) 
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Figure 0.5. (a) The Frozen Wave with Q = 0.999ojo /c and N = 45, approximately reproducing 
the chosen longitudinal pattern represented by eq. d0.77t , (b) A different Frozen wave, now with 
Q — 0.995wo /c (but still with N = 45) forwarding the same longitudinal pattern. We can observe 
that in this case (with a lower value for Q) a higher transverse localization is obtained. 



with transverse size, or width, Ap rn w 2A/k pm . On the basis of the expected convergence 
of series (10.741 1. we can estimate the width of the transverse spot of the resulting beam as 
being 

2.4 2.4 

A <° w Z = / 2 i 1 m ' (0 - 78) 

k pm=0 ^ojI/c 2 -Q 2 

which is the same value as that for the transverse spot of the Bessel beam with m = 
in superposition (10.741 . Relation (10.78b can be useful: Once we have chosen the desired 
longitudinal intensity pattern, we can choose even the size of the transverse spot, and use 
relation HO. 78\) for evaluating the needed, corresponding value of parameter Q. 

For a more detailed analyzis concerning the spatial resolution and residual intensity of 
the Frozen Waves, we refer the reader to ref.[17]. 



0. 5. 1. 1 Increasing the control on the transverse shape by using higher-order 
Bessel beams 

Here, we are going to argue that it is possible to increase even more our control on 
the transverse shape by using higher-order Bessel beams in our fundamental superposition 
(10741 

This new approach can be understood and accepted on the basis of simple and intuitive 
arguments, which are not presented here, but can be found in ref.[17]. A brief description 
of that approach follows below. 

The basic idea is obtaining the desired longitudinal intensity pattern not along the axis 
p = 0, but on a cylindrical surface corresponding to p = p 1 > 0. 

To do this, we first proceed as before: Once we have chosen the desired longitudinal inten- 
sity pattern F(z), within the interval < z < L, we calculate the coefficients A m as before, 

1. e., A m = (1/L) J Q L F(z) exp(—i2Ttmz/L) dz, and k pm = d uj 2 — (Q + 2Trm/L) 2 . 
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Afterwards, we just replace the zero-order Bessel beams Jo(k pm p), in superposition 
( 10.74b . with higher-order Bessel beams, Jfj,(k pm p), to get 

N 

*(p,z,t) = e -^ ot e tQz AmMkpmP)e i3 * m ' , (0.79) 

m=-N 

With this, and based on intuitive arguments [17], we can expect that the desired longitu- 
dinal intensity pattern, initially constructed for p — 0, will approximately shift to p = p' , 
where p 1 represents the position of the first maximum of the Bessel function, i.e the first 
positive root of the equation (d J/j,(kp m —o p) /dp)\ p > = 0. 

By such a procedure, one can obtain very interesting stationary configurations of field 
intensity, as "donuts", cylindrical surfaces, and much more. 

In the following example, we show how to obtain, e.g., a cylindrical surface of stationary 
light. To get it, within the interval < z < L, let us first select the longitudinal intensity 
pattern given by eq. (10.771) . with U = L/2-Azwdl 2 = L/2 + Az, and with Az = L/300. 
Moreover, let us choose L = 0.05 m, Q = 0.998 ujq/c, and use N = 150. 

Then, after calculating the coefficients A m by eq. (10.731) . we have recourse to su- 
perposition (10.79b . In this case, we choose p, = 4. According to the previous dis- 
cussion, one can expect the desired longitudinal intensity pattern to appear shifted to 
p 1 ps 5.318/fcp m= o = 8.47 pm, where 5.318 is the value of k pm -o p for which the Bessel 
function J4 (k p m= o p) assumes its maximum value, with k p TO= o = yo^ — Q 2 - The figure 
I0.6.l below shows the resulting intensity field. 

In Fig l0.6.f a) the transverse section of the resulting beam for z — L/2 is shown. The 
transverse peak intensity is located at p = 7.75 pm, with a 8.5% difference w.r.t. the 
predicted value of 8.47 pm. Figure lOlxT b) shows the orthogonal projection of the resulting 
field, which corresponds to nothing but a cylindrical surface of stationary light (or other 
fields). 





Figure 0.6. (a) Transverse section at z = L/2 of the considered higher-order FW. (b) Orthogonal 
projection of the three-dimensional intensity pattern of the same higher-order FW. 
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We can see that the desired longitudinal intensity pattern has been approximately ob- 
tained, but, as wished, shifted from p = to p = 7.75 pm, and the resulting field resembles 
a cylindrical surface of stationary light with radius 7.75 pm and length 238 pm. Donut-like 
configurations of light (or sound) are also possible. 

0.5.2 Stationary wavefields with arbitrary longitudinal shape in absorbing 
media: Extending the method. 

When propagating in a non-absorbing medium, the so-called nondiffracting waves main- 
tain their spatial shape for long distances. However, the situation is not the same when deal- 
ing with absorbing media. In such cases, both the ordinary and the nondiffracting beams 
(and pulses) will suffer the same effect: an exponential attenuation along the propagation 
axis. 

Here, we are going to make an extension[10] of the method given above to show that, 
through suitable superpositions of equal-frequency Bessel beams, it is possible to obtain 
nondiffracting beams in absorbing media, whose longitudinal intensity pattern can assume 
any desired shape within a chosen interval < z < L of the propagation axis z. 

As a particular example, we obtain new nondiffracting beams capable to resist the loss 
effects, maintaining amplitude and spot size of their central core for long distances. 

It is important to stress that in this new method there is no active participation of the 
material medium. Actually, the energy absorption by the medium continues to occur nor- 
mally, the difference being in the fact that these new beams have an initial transverse field 
distribution, such to be able to reconstruct (even in the presence of absorption) their central 
cores for distances considerably longer than the penetration depths of ordinary (nondiffract- 
ing or diffracting) beams. In this sense, the present method can be regarded as extending, 
for absorbing media, the self-reconstruction properties [27] that usual Localized Waves are 
known to possess in loss-less media. 

In the same way as for lossless media, we construct a Bessel beam with angular frequency 
uj and axicon angle 9 in the absorbing materials by superposing plane waves, with the same 
angular frequency u, and whose wave vectors lie on the surface of a cone with vertex angle 
9. The refractive index of the medium can be written as n(u)) = tih^uj) + inj(ui), quantity 
nji being the real part of the complex refraction index and n/ the imaginary one, responsible 
for the absorbtion effects. With a plane wave, the penetration depth S for the frequency u 
is given by S = l/a = c/2umi, where a is the absorption coefficient. 

In this way, a zero-order Bessel beam in dissipative media can be written as ip = 
Jo(k p p)exp(if3z)exp(—icLit) with (3 = n(oj)uj cos 9 / c = njiojcos9/c + iniujcos9/c = 
Pr + ifli', k p = njiujsm9/c + iniu;sm9/c = k p R + ik p j, and so k 2 p = n 2 uj 2 /c 2 — 1 . 
In this way,] it results ip = Jo((k p R + ik p ])p)exp(iPnz)exp(—iujt)exp(~Pjz), where 
/3r, k p R are the real parts of the longitudinal and transverse wave numbers, and /?/, k p j 
are the imaginary ones, while the absorption coefficient of a Bessel beam with axicon 
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angle 9 is given by ag — 2/3/ = 2n/wcos0/c, its penetration depth being Sg = 1/ag — 
c/2wn/ cos 9. 

Due to the fact that k p is complex, the amplitude of the Bessel function Jo(k p p) starts 
decreasing from p = till the transverse distance p = l/2k p ], and afterwards it starts 
growing exponentially. This behavior is not physically acceptable, but one must remember 
that it occurs only because of the fact that an ideal Bessel beam needs an infinite aperture 
to be generated. However, in any real situation, when a Bessel beam is generated by finite 
apertures, that exponential growth in the transverse direction, starting after p = l/2k p j, 
will not occur indefinitely, stopping at a given value of p. Let us moreover emphasize that, 
when generated by a finite aperture of radius R, the truncated Bessel beam[17] possesses a 
depth of field Z = R/ tan 9, and can be approximately described by the solution given in 
the previous paragraph, for p < R and z < Z. 

Experimentally, to guarantee that the mentioned exponential growth in the transverse 
direction does not even start, so as to meet only a decreasing transverse intensity, the radius 
R of the aperture used for generating the Bessel beam should be R < l/2fc p /. However, as 
noted by Durnin et al., the same aperture has to satisfy also the relation R > 2w/k pR . From 
these two conditions, we can infer that, in an absorbing medium, a Bessel beam with just 
a decreasing transverse intensity can be generated only when the absorption coefficient is 
a < 2/A, i.e., if the penetration depth is 6 > A/2. The method developed in this paper does 
refer to these cases, i.e. we can always choose a suitable finite aperture size in such a way 
that the truncated versions of all solutions presented in this work, including the general one 
given by eq.(6), will not present any unphysical behavior. Let us now outline our method. 

Consider an absorbing medium with the complex refraction index n(uj) = n R (oj) + 
ini(uj), and the following superposition of 2N + 1 Bessel beams with the same frequency 
ur. 



n 

9(p,Z,t) = A mJo((k pRm +ik pIm )p) e i0n m z e -iut e -f} lmZ ^ (Q g()) 

m=-N 

where the m are integer numbers, the A m are constant (yet unknown) coefficients, quantities 
/?R m and k p R m (/3/ m and k p j m ) are the real (the imaginary) parts of the complex longitudinal 
and transverse wave numbers of the m-th Bessel beam in superposition (10.801) : the following 
relations being satisfied 

k 2 pm = n 2 ^-p 2 m (0.81) 

^ = I* (0.82) 

where (3 m = (3 Rm + i(3 Im , k Pm = k pRm + ik pIm , with k pRm /k pIm = n R /m. 

Our goal is now to find out the values of the longitudinal wave numbers f3 m and the 
coefficients A m in order to reproduce approximately, inside the interval < z < L (on the 
axis p = 0), a freely chosen longitudinal intensity pattern that we call |F(z)| 2 . 
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The problem for the particular case of lossless media[16,17], i.e., when m — — > 
[3j m — 0, was solved in the previous subsection. For those cases, it was shown that the 
choice — Q + 2itm/ L, with A m = F(z)exp(—i27rmz/L)/L dz can be used to 
provide approximately the desired longitudinal intensity pattern \ F(z) | 2 on the propagation 
axis, within the interval < z < L, and, at the same time, to regulate the spot size of the 
resulting beam by means of the parameter Q, which parameter can be also used to obtain 
large field depths and also to inforce the linear polarization approximation to the electric 
field for the TE electromagnetic wave (see details in refs.[16,17]). 

However, when dealing with absorbing media, the procedure described in the last para- 
graph does not work, due to the presence of the functions exp(— /3j m z) in the superposition 
( 10.80I I. since in this case that series does not became a Fourier series when p = 0. 

On attempting to overcome this limitation, let us write the real part of the longitudinal 
wave number, in superposition ( 10.801 ). as 



with 



P Rm = Q + ^ (0.83) 



0<Q + ^<n^ (0.84) 
L c 

where this inequality guarantees forward propagation only, with no evanescent waves. 
In this way the superposition d0.801 > can be written 

N 

*{p,z,t) = e- i " t e i Q z J2 A mJo((k pRm +ik pIm )p) j^* e -Pr m * f (0.85) 

m=-N 

where, by using (10.82b . we have /3j m = (Q + l-nrnj L)nj /ur, and k Pm = k pRm + ik p i m 
is given by (10.81b . Obviously, the discrete superposition (10.85b could be written as a 
continuous one (i.e., as an integral over /3_R m ) by taking L — > oo, but we prefer the discrete 
sum due to the difficulty of obtaining closed-form solutions to the integral form. 

Now, let us examine the imaginary part of the longitudinal wave numbers. The minimum 
and maximum values among the f3j m are (/3/) m i n = (Q — 2irN j L)ni/nn and (/3/) ma x = 
(Q + 2TrN/L)ni/nR, the central one being given by (3j = (j3i) m =o — Qni/riR. With 
this in mind, let us evaluate the ratio A = [(/?/ ) ma x — (Pi)min]/ Pi — 4:TtN/LQ. 

Thus, when A << 1, there are no considerable differences among the various /3j m , 
since /3/ m ~ (3j holds for all m. In the same way, there are no considerable differences 
among the exponential attenuation factors, since exp(— @i m z) ~ exp(— fijz). So, when 
p = the series in the r.h.s. of eq. (10.85l ) can be approximately considered a truncated 
Fourier series multiplied by the function exp(— Pjz) and, therefore, superposition ( 10.851 ) 
can be used to reproduce approximately the desired longitudinal intensity pattern |i^(z)| 2 
(on p = 0), within < z < L, when the coefficients A m are given by 

A m = ~ f F(z)e^' z e- t2 ^ z dz, (0.86) 
L Jo 
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being necessary the presence of the of the factor exp(/3 7 z) in the integrand to compensate 
for the factors exp(— (3i m z) in superposition ( 10.85b . 

Since we are adding together zero-order Bessel functions, we can expect a good field 
concentration around p = 0. 

In short, we have shown in this Section how one can get, in an absorbing medium, a sta- 
tionary wave-field with a good transverse concentration, and whose longitudinal intensity 
pattern (on p = 0) can approximately assume any desired shape \F(z)\ 2 within the pre- 
determined interval < z < L. The method is a generalization of a previous one[16,17] 
and consists in the superposition of Bessel beams in EO.( l0.85b . the real and imaginary 
parts of their longitudinal wave numbers being given by eqs.( 10.83t and ( 10.821 ). while their 
complex transverse wave numbers are given by eq.( IQ.81l ). and, finally, the coefficients of 
the superposition are given by eq. ( 10.86b . The method is justified, since AttN/LQ << 1; 
happily enough, this condition is satisfied in a great number of situations. 

Regarding the generation of these new beams, given an apparatus capable of generat- 
ing a single Bessel beam, we can use an array of such apparatuses to generate a sum of 
them, with the appropriate longitudinal wave numbers and amplitudes/phases [as required 
by the method], thus producing the desired beam. For instance, we can use [16, 17] a laser 
illuminating an array of concentric annular apertures (located at the focus of a convergent 
lens) with the appropriate radii and transfer functions, able to yield both the correct longi- 
tudinal wave numbers (once a value for Q has been chosen) and the coefficients A n of the 
fundamental superposition (10.85b . 

0.5.2. 1 Some Examples For generality's sake, let us consider a hypothetical medium 
in which a typical XeClexcimer laser (A = 308nm — > u> = 6.12xl0 15 Hz)hasapenetration 
depth of 5 cm; i.e. an absorption coefficient a = 20m _1 , and therefore rii — 0.49 x 10~ 6 . 
Besides this, let us suppose that the real part of the refraction index for this wavelength is 
riR = 1.5 and therefore n — ur + ini = 1.5 + i0.49 x 10~ 6 . Note that the value of 
the real part of the refractive index is not so important for us, since we are dealing with 
monochromatic wave fields. 

A Bessel beam with ui = 6.12 x 10 15 Hz and with an axicon angle 9 = 0.0141 rad 
(so, with a transverse spot of radius 8.4 pm), when generated by an aperture, say, of radius 
R = 3.5 mm, can propagate in vacuum a distance equal to Z = Rj tan^ = 25 cm while 
resisting the diffraction effects. However, in the material medium considered here, the 
penetration depth of this Bessel beam would be only z p = 5 cm. Now, let us set forth two 
interesting applications of the method. 

First Example: Almost Undistorted Beams in Absorbing Media. 

Now, we can use the extended method to obtain, in the same medium and for the same 
wavelength, an almost undistorted beam capable of preserving its spot size and the intensity 
of its central core for a distance many times larger than the typical penetration depth of an 
ordinary beam (nondiffracting or not). 
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With this purpose, let us suppose that, for this material medium, we wish a beam (with 
to = 6.12 x 10 15 Hz) that maintains amplitude and spot size of its central core for a distance 
of 25 cm, i.e. a distance 5 times greater than the penetration depth of an ordinary beam 
with the same frequency. We can model this beam by choosing the desired longitudinal 
intensity pattern |F(z)| 2 (on p = 0), within < z < L, 

F(z) = { 1 im °< Z < Z (087) 
[0 elsewhere, ' 

and by putting Z = 25 cm, with, for example, L = 33 cm. 

Now, the Bessel beam superposition ( 10.851 ) can be used to reproduce approximately this 

intensity pattern, and to this purpose let us choose Q = 0.9999o;/c for the f3n m in ( 10.83b . 

and N — 20 (note that, according to inequality ( 10.841 ). TV could assume a maximum value 

of 158.) 

Once we have chosen the values of Q, L and N, the values of the complex longitudinal and 
transverse Bessel beams wave numbers happen to be defined by relations (10.831 ). (10.821 ) and 
(10.811) . Eventually, we can use eq. (10.86l) and find out the coefficients A m of the fundamental 
superposition (10.851) . that defines the resulting stationary wave-field. 

Let us just note that the condition AttN/LQ << 1 is perfectly satisfied in this case. 

In Fig. (10. 7. 1) ("a) we can see the 3D field-intensity of the resulting beam. One can see 
that the field possesses a good transverse localization (with a spot size smaller than 10 /im), 
it being capable of maintaining spot size and intensity of its central core till the desired 
distance (a better result could be reached by using a higher value of N). 




z (mm) 



(a) (b) 

Figure 0.7. (a) Three-dimensional field-intensity of the resulting beam, (b) The resulting beam, in 
an orthogonal projection and in logaritmic scale. 

It is interesting to note that at this distance (25 cm), an ordinary beam would have got 
its initial field-intensity attenuated 148 times. 

As we have said in the Introduction, the energy absorption by the medium continues to 
occur normally; the difference is that these new beams have an initial transverse field dis- 
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tribution sophisticated enough to be able to reconstruct (even in the presence of absorption) 
their central cores till a certain distance. For a better visualization of this field-intensity 
distribution and of the energy flux, Fig. (10.7.K b) shows the resulting beam, in an orthogo- 
nal projection and in logarithmic scale. It is clear that the energy comes from the lateral 
regions, in order to reconstruct the central core of the beam. On the plane z = 0, within 
the region p < R = 3.5 mm, there is a uncommon field intensity distribution, being very 
dispersed instead of concentrated. This uncommon initial field intensity distribution is re- 
sponsible for constructing the central core of the resulting beam and for its reconstruction 
till the distance z = 25 cm. Due to absorption, the beam (total) energy flowing through 
different z planes, is not constant, but the energy flowing IN the beam spot area and the 
beam spot size itself are conserved till (in this case) the distance z = 25 cm. 

Second Example: Beams in absorbing media with a growing longitudinal field intensity. 

Considering again the previous hypothetical medium, in which an ordinary Bessel beam 
with 9 = 0.0141 rad and uj = 6.12 x 10 15 Hz has a penetration depth of 5 cm, we aim to 
construct now a beam that, instead of possessing a constant core-intensity till the position 
z = 25 cm, presents on the contrary a (moderate) exponential growth of its intensity, till 
that distance (z — 25 cm). 

Let us assume we wish to get the longitudinal intensity pattern |F(z)| 2 , in the interval 
< z < L, 



with Z = 25 cm and L = 33 cm. 

Using again Q = 0.9999w/c, N = 20, we CAN proceed as in the first example, 
calculating the complex longitudinal and transverse Bessel beams wave numbers and finally 
the coefficients A m of the fundamental superposition (I0.851 l. 

In Fig. (10.8.1 1 we can see the 3D field-intensity of the resulting beam. One can see that the 
field presents the desired longitudinal intensity pattern with a good transverse localization 
(a spot size smaller than 10 /im). 

Obviously, the amount of energy necessary to construct these new beams is greater than 
that necessary to generate an ordinary beam in a non-absorbing medium. And it is also clear 
that there is a limitation on the depth of field of these new beams. In the first example, for 
distances longer than 10 times the penetration depth of an ordinary beam, besides a greater 
energy demand, we meet the fact that the field-intensity in the lateral regions would be even 
higher than that of the core, and the field would loose the usual characteristics of a beam 
(transverse field concentration) 




exp(>/Z) for < z < Z 
elsewhere, 



(0.88) 
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Figure 0.8. Three-dimensional field-intensity of the resulting beam, in an absorbing medium, with 
a growing longitudinal field intensity. 
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